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ABSTRACT 

In this paper a multi item EOQ model with stock dependent demand for deteriorating items is considered in fuzzy 
environment. Different inventory costs and the amount of investment are represented as triangular fuzzy numbers The model 
has been solved by Robust ranking method, fuzzy optimization technique (FOT) and intuitionistic fuzzy optimization 
technique (IFOT). Nearest interval approximation method is used in FOT and IFOT. Pareto optimality test is done for fuzzy 
optimization and intuitionistic fuzzy optimization techniques. The methods are illustrated with a numerical example. 
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I. INTRODUCTION 

Inventory problems are common in manufacturing, maintenance service and business operations in general. 
Often uncertainties may be associated with demand, various relevant costs. Generally demand rate is considered to be 
constant, time dependent, ramp type and selling price dependent. However in present competitive market stock- 
dependent demand plays an important role in increase its demand. Deterioration is one of the important factors in 
inventory system. Some items like food grains, vegetables, milk, eggs etc. deteriorating during their storage time and 
retailer suffers loss. In an inventory system available storage space, budget, number of orders etc. are always limited 
hence multi-item classical inventory models under these constraints have great importance. 

In conventional inventory models, uncertainties are treated as randomness and are handled by probability 
theory. However in certain situations, uncertainties are due to fuzziness and in such cases the fuzzy set theory 
introduced by Zadeh [18] is applicable. Many researches introduced the concept of ranking function [] for fuzzy 
parameter which gives a representative value for it. Using ranking function fuzzy parameters can be easily changed 
into crisp one. 

The fuzzy parameters changing into appropriate interval numbers are introduced by Gregorzwshi [5],Susovan 
chakrabortly et.al [15,16] proposed a method to solve an EOQ model using the concept of interval arithmetic. In fuzzy 
optimization the degree of acceptance of objectives and constraints are considered here. Now a days different 
modification and generalized form of fuzzy set theory have appeared intuitionistic fuzzy set (IFS). The concept of an 
IFS can be viewed as an alternative approach to define a fuzzy set in case where available information is not sufficient 
for the definition of an imprecise concept by means of a conventional fuzzy sets. 
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In this paper, a multi-item inventory model of deteriorating item with stock-dependent demand is formulated in 
crisp and fuzzy environment. Here the objective is to find an EOQ to maximize the profit. The problem is solved by 

• Robust’s ranking technique 

• Fuzzy optimization technique 

• Intuitionistic fuzzy optimization technique 

The model is illustrated numerically and the results are obtained from different methods. 

II.PRELIMINARIES 


Definition (1): Atanassov [20], Let X={xl x 2 , x n } be a finite universal set. An Atanassov’s intuitionistic fuzzy 

set (IFS) in a given universal set X is an expression A is given by £ _ [i „ (x .) v ■ (x )\ x & X } 

Where the functions ^:X->[0,1] i-e. x. . e A — > fi . (x. ) :e [0, 1] and T A :X->[0,11 i.e. 
X. € X — > V- ( x i ) € 1 0, 1 1 define the degree of membership and the degree of non-membership respectively of an 
element X t E X satisfy the condition: for every x . <= x , o < //-(jc) + v- (jc)< 1 - 


Definition (2): Let A and B be two Atanassov’s IFSs in the set X. The intersection of A and B is defined as 


follows: 


An B = min (ju A (x ; ), ju - (x,.)), max (v A (x , ), v. {x t )))|x ; e X }■ 


Definition (3): Let 9i be the set of all real numbers. An interval, moore[9], may be expressed as 

a =[a L ,a R ] ={x:a L <x<a R ,a L e%a R e 91 } 


( 1 ) 


Where a L and a R are called the lower and upper limits of the interval a , respectively. 

If a L = a R then a = [a L , a R ] is reduced to a real number a, where a= a L = a R . alternatively an interval a can be 
expressed in mean-width or center-radius form as Cl = (m(a) , .where m(a ) = — [a L + a R ) and 

w(a) = — (ci R — a L ) are respectively the mid-point and half-width of the interval a .The set of all interval 

numbers in Ji is denoted by 1(91). 

Definition (4): (Triangular Fuzzy Number) 

For a triangular fuzzy number A(x), it can be represented by A(a,b,c:l) with membership function Pa(x) given by 


M a ( x ) = 


( x - a) 


a < x < b 


(b - a) 

1 \ x = b 

(c ~ x) 

(c - b) 

0 ; otherwise 


; b < x < c 
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Definition (5): (a-cut of a Fuzzy Number) 

The a-cut of a fuzzy number A(x) is defined as A(a)={x:p(x)>a, ae[0,l] } 

Definition (6): Nearest Interval Approximation 


According to Gregorzewski [5] we determine the interval approximation of a fuzzy number as: Let A 
=(a b a 2 ,a 3 ) be an arbitrary triangular fuzzy number with a a-cut [A L (a),A R (a)] 

Then by nearest interval approximation method, the lower limit Cl and upper limit C R of the interval are 

l l 


C L = J A L {a)da = J[a 1 + (a 2 -ajajcfa = 


a { + a,, 


t t 

C R =^A R (a)da = ^\ji 3 -(a 3 -a 2 )a~ 1 \fla = 


^2 I 


( 2 ) 


I ^2 $2 I $3 


Therefore, the interval number considering A as triangular fuzzy number is 

Definition)?) : Pareto-Optimality 

x e X is said to be a Pareto optimal solution if there does not exist another XG X S uch that /* ( x ) — fk ( x ) 
forallP=l,2, p and f j( x ) < 'fj( x ) for at least one p } 

III. ASSUMPTION 


• The scheduling period is constant and no lead time. 

• Replenishment rate is infinite. 

• Selling price is known and constant. 

• Demand rate is stock dependent. 

• Shortages are not allowed. 

• Deteriorating rate is age specific failure rate. 

IV.NOTATIONS 

Tj : Time period for each cycle for the i th item. 

Ri : Demand rate per unit time of i th item. | R = a.+h.q,] 

0i : Deteriorating rate of i th item. 

Qi(t) : inventory level at time t of i th item. 

C H : Total Holding cost. 

Cn :Holding cost per unit of i th item. 
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C 3i : Setup cost for i th item. 

S d ; : Total deteriorating units of i th item. 

P; : Selling price per unit of i th item. 

Qi :Initial stock level of i th item. 

PF(Qi) : Total profit of i th item. 

N :Number of items. 

(wavy bar (~) represents the fuzzification of the parameters ) 

V. MATHEMATICAL FORMULATION 

5.1 Crisp Model 


As Qi(t) is the inventory level at time t of i th item, then the differential equation describing the state of inventory is 


given by 


— Q i {t) + 6 i Q i {t) = -{a i +b i Q i {t)) 
dt 


0<t <T t 

Solving the above differential equation using boundary condition Qi(t)= Q, at t=0, we get 


2 ,( 0 =- 


(0, + b,) 


Q, 


(8, + b,) 




(3) 


And using boundary condition Qi(t)= 0 at t=T, 


••• T = log 1 1 + 

(0i + b i) \ a , 


(4) 


The holding cost of i th item in each cycle is 
Ch = Ch Gj (Qi) 

Where, 


(5) 


o,(e,)= J 


< hdq t 


a, + (0, + b j )q l 


Q, 


- + ■ 


-log( 1 + 


(O.+b.lQ, 


(0,+b,) (e+bf 

By neglecting the higher power terms, we get 


G i {Q l ) = ^\\-2 {9i + bi)Qi 
2a. 3a, 


The total number of deteriorating units of the i th item is 
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The net revenue for the i th item is 


N{Q i ) = {P i -C i )Q i -P i .S di {Q i ) 
N{Q i ) = {P i -C i )Q i -P i e i .G i {Q i ) 

The profit of i th item is 


PF(Q l ) = N(Q i )-C u G i (Q i )~C 


3 !’ 


i= 1 ,2, n. 


pm )=(p- c, )q, - pa (Qi ) - c i; g, (a ) - c 3i 


i= 1 ,2, n. 


(7) 


PF(Qi ) = (P- q )Q - (C u + PA )Gj (a ) ■ - C 3i , i=L,2, n. 


Hence the problem is 

Max = X[(^- C «)a-(C u + W)G,(e,)-C J ,] (8) 

(=1 


Qi > 0, i=l,2, n. 

5.2. Fuzzy Model 

When above profit, cost, selling price, deteriorating rate and total budget become fuzzy, the said crisp model is 
transform to 


M ax PF = £ [(/> - C,)Q, - (C„ + PfiJG.lQ,) - C 3t ] W 

j=i 

Qi>0, i= 1 ,2, n 

VI.MATHEMATICAL ANALYSIS 

• Robust Ranking Technique 

Given a convex fuzzy number a , the Robust’s ranking index is defined by 
1 

R ( d ) = j0.5(a a L ,a a u )da 

0 

5 

Where {^Cl a L , Cl a U ) is the a-level cut of the fuzzy number b, . 


In this paper we use this method of ranking for the fuzzy numbers. The Robust’s ranking index R(£I) gives the 
representative value of the fuzzy number Cl . Hence using this we can defuzzify into crisp one. 

Robust’s ranking technique [10,14] which satisfies compensation, linearity and additivity properties and provides 
results which are consistent with human intuition. 

We apply Robust’s ranking method[10,14] to defuzzify the fuzzy objective. 
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Taking P ; i= 1 ,2, n as fixed in (9) and Using this ranking function the fuzzy model (9) becomes, 

M ax PF* = ^ [(/* - R{C i ))Q j - R(C u )G l (Q l ) - P i R(0 i )G i (Q i ) - /?(C 3i )] 

i = 1 L 

• B. Fuzzy optimization Technique 

Using nearest interval approximation [5] to the objective and Robust’s ranking [10,14] to the constant. The 
proposed model can be stated as 

max imize { PF L (Q i ), PF R (Q i )} (10) 

Here the interval valued problem (10) is represented as 

max imize {PF l ( Q . ), PF C (Q, )PF K (Q i )} PF C =^[PF,+ PF R ] (11) 

To solve the above multi-objective problem (11) we have used the following fuzzy programming technique. 

Step 1: Solve the multi-objective profit function as a single objective Profit function using one objective at a time 
and ignoring all others. 

Step 2: From the results of step 1, determine the corresponding values for every objective at each solution 

derived. 

Step 3: From step 2, obtain the lower bounds and upper bounds for each objective functions. Set U k =max (PF k ) 
and L k = min iPF k ) K=L,C,R. 

Step 4: Define fuzzy linear membership function ^ FpF k \ k = L,C, R^J for each objective function is defined by 

0 ; PF k < L k 
PF -L 

MrF k = \ \ - ; h ^ PF k £ U , ( 12 ) 

d k 

1 ; PF k >U k 

Step 5: After determining the linear membership function defined in( 12) for each objective functions following 
the problem (11) can be formulated an equivalent crisp model 

Max a, 

a < ju PFk (Q);k = L,C,R-, O<a<l,0,. >0 

Step 6: Now the above problem can be solved by a non-linear optimization technique and optimal solution of a 
(say a*) is obtained. 

Step 7: Now after obtaining a*, the decision maker selects the objective function, from among the objective 
functions PF L , PF C , PF R . If he selects PF C 
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Max PF C , 


PF l ^ m L , PF C < m c , PF r < m R 


(13) 


Where m L =U L -a*d L ;m c = U c -cx*d c ;m R = U R - a * d R ; 

Step8: Pareto-optimal solution 

Now after deriving the optimum decision variables, Pareto-optimality test is performed according to [16], let the 
decision vectors Q t and the optimum values PF L = PF L (Q ), PF C = PF C (Q ), PF R = PF R (Q ), are obtained 
from (13). With these values, the following problem is solved using a non-linear optimization technique. 


min V = (e L +e c +e R ) 


Subject 


(14) 


to: 


PF +<= = PF PF +e =PF PF -l-e = PF 

1 1 L' L 1 1 L ’ 1 1 C' C 11 C’ 11 R^^R 1 1 R ’ 


e L ,e c ,e R >0,Q>0,a> j3 
a+ /3<l,a,/3>0 


__ ^ssj: __ w 

, PF l , PF C , PF r are called strong Pareto Optimal solution provided V 
is very small otherwise it is called weak Pareto solution. 

• C. IF Programming Technique 

To solve multi-objective maximization problem given by (14), we have used the following IF programming 
technique. 


For each of the objective functions P F L (Q), P F c ( Q), PF ;; ( Q ) , we first find the lower bounds L L , L c , L R 
(best values) and the upper bounds JJ L ,U c ,U R (worst values), where L L , L c , L R are the aspired level achievement and 
U L , U c , U R are the highest acceptable level achievement for the objectives p F L (Q ), P F c (Q), PF ;( (Q ) respectively 

and d k = U k — L k is the degradation allowance for objective PF k (Q ) , k=L,C,R. Once the aspiration levels and 

degradation allowance for each of the objective function has been specified, we formed a fuzzy model and then transform 
the fuzzy model into a crisp model. The steps of intuitionistic programming technique is given below. 

Step 1: Solve the multi-objective profit function a single objective Profit function using one objective at a time 
and ignoring all others. 

Step 2: From the results of step 1, determine the corresponding values for every objective at each solution 

derived. 


Step 3: From step 2, obtain the lower bounds and upper bounds for each objective functions. Set Uk=max (PF k ) 
and L k = min (PFk) , K=L,C,R. 


Step 4: Define membership function 



L,C,R j , and non membership function ( P /v _- ; k = L,C,R 
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for each objective function. A linear membership function is defined by 


H-pf. 


; PF K <L k 


P F k ~L k 
d. 


■L k <PF k <U k 


■,pf k >u k 


A Linear non- membership function is defined by 


0 

1- 

1 


■PF K <L t 


r PF iZ L t '' 
V d k J 


,L k <PF k <U k 


■,pf k >u k 


(15) 


Step 5: After determining the membership and non-membership function defined in(14) for each objective 
functions following the problem (13) can be formulated an equivalent crisp model on the basis of definition 2 of this paper 
as 

Max a . min P 

a<jU PFt (x) ;k = L,C,P 

j3>v PFk (x ) ;k = L,C,R 

a>/3. a + J3 <\;a, J3 > 0 


0 < a < 1 , Q t > 0 

Step 6: Now the above problem can be solved by a non-linear optimization technique and optimal solution of a 
(say a*) and p, (say P*) are obtained. 

Step 7: Now after obtaining a* and P*. the decision maker selects the objective function, from among the 
objective functions PF L , PF C , PF R . If he selects PF C 

MaxPF c 

PF l < m L , PF C < m c , PF R < m R 
p F L >n L ,PF c >n c ,PF R >n R > 

Q > 0,a > fd\ 
a+jB<l; a,jB> 0 

m L = L l + ad L ; n L =U L — j3d L 
Where m c = L ( + ad c ;n c = U C — J3d c 
m R = L r + ad R ; n R = U R — j8d R 
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Step 8: Pareto-optimal solution 

Now after deriving the optimum decision variables, Pareto-optimality test is performed according to [16], let the 
decision vector Q and the optimum values PF L = PF L (Q ), PF ( = PF C (Q ), PF R = PF R (Q ) , are obtained 
from (14). With these values, the following problem is solving using a non-linear optimization technique. 


Subject 


(16) 


minV = (e L +e c +e s ) 

PF +e = PF* PF +€ =PF * PF +e =PF* 

to: > 

e L ,e c ,e R >0,Q>0,a> j8 
a+ fi <l;a,j6>0 


The optimal solution of (16), say {J , P F, , PF C , PF K are called strong Pareto Optimal solution provided V 
is very small otherwise it is called weak Pareto solution. 

VII.NUMERICAL EXAMPLE 


For all models let us assume 
For n=2, 

p 1= 10, C|=7, a 1= l 10, b |=0.5, fl^O.025, P 2 = 10, C 2 =6.75, a 2 =100, b 2 =0.5, 9 2 =0.03, C„=2, C 12 =2.2, C 3 i=65, 
C 32 =50, B=1800. 

• Taking C n =(1.95,2,2.05); C 12 =(2.15,2.20,2.25); C 3] =(55,65,75); C 32 =(40,50,60); 6 X 

=(0.02,0.025,0.03) 0 2 =(0.025,0.03,0.035); C x =(6.5, 1, 7.5); C 2 =(6.25,6.75,7.25); B 

=(1700,1800,1900) 

max PF* = 3g, + 3.25g 2 - 0.0102262i 2 - 0.0125e 2 2 - 1 

The optimum results are 


Table 1 


Max PF* 

Qi* 

q 2 * 

316.27738 

146.68492 

130 


• Considering all values as same as case (A) and also taking P ] =(9.5,10,1 1); P n =(9.75, 10,1 1);. 


Table 2 



PF, 

PF r 

pf r 

PF l 

184.67881 

162.78795 

90.677814 

PF C 

342.71953 

363.87741 

339.55730 

pf r 

500.75962 

564.96589 

588.43534 


The optimum value of a is 0.75002106 
Optimum results when PF L is chosen as the most important objective function 
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Table 3 


pf l * 

PF C * 

pf r * 

Qi* 

q 2 * 

161.18053 

357.79789 

554.4140663 

132.77717 

152.94231 


Pareto Optimal Solution 


Table 4 


V 

pp L ** 

PF C ** 

PF r ** 

Ox** 

q 2 ** 

0.70661E-05 

161.18052 

357.79788 

554.41407 

132.77717 

152.94232 


• Considering all values as same as case (B) and also taking P { =(9.5,10,1 1); P 2 =(9.75, 10,1 1);. 


The optimum value of a and P 


Table 5 


a 

P 

0.75002106 

0.24997894 


Optimum results when PF L is chosen as the most important objective function 

Table 6 


PF, * 

PFr* 

PFr* 

Qi* 

02* 

161.18054 

363.84990 

566.5182 

153.74993 

140.95712 


Pareto Optimal Solution 


Table 7 


V 

pf l ** 

PF r ** 

PF r ** 

Ox** 

q 2 ** 

0.33021E-04 

161.18053 

363.84986 

566.51819 

153.71317 

140.99336 


VIII.CONCLUSIONS 

In this paper, a real life inventory problem under the investment constraint environment has been proposed and 
solved by three different methods namely (i) Robust’s ranking (ii) fuzzy optimization (FOT) (iii) intuitionistic fuzzy 
optimization (IFOT). Different inventory costs and the investment amount are considered as triangular fuzzy numbers 
nearest interval approximation method is used in FOT and IFOT. Linear membership function is taken for fuzzy objective. 
Robust’s ranking method gives the less profit amount when compared with fuzzy optimization and intuitionistic fuzzy 
optimization. FOT and IFOT yields better results. Also strong Pareto optimality is obtained for FOT and IFOT. 
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